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We investigate collective oscillations of non-degenerate clouds of 87 Rb atoms as a function of 
density in an elongated magnetic trap. For the low-lying M = monopole-quadrupole shape os- 
cillation we measure the oscillation frequencies and damping rates. At the highest densities the 
mean-free-path is smaller than the axial dimension of the sample, which corresponds to collisionally 
hydrodynamic conditions. This allows us to cover the cross-over from the collisionless to the hy- 
drodynamic regime. The experimental results show good agreement with theory. We also analyze 
the influence of trap anharmonicities on the oscillations in relation to observed temperature depen- 
dencies of the dipole and quadrupole oscillation frequencies. We present convenient expressions to 
quantify these effects. 
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I. INTRODUCTION 

Collisional hydrodynamics has gradually become an 
important issue for the understanding of experiments 
with dilute quantum gases. When the atomic mean- free- 
path is smaller than the characteristic dimensions of typ- 
ical elongated atomic clouds, the gas properties depend 
on the local density field and exhibit collisional hydrody- 
namics rather than the collisionless dynamics of a nearly 
ideal gas 0, Q- For Bose gases and Bose- Fermi mix- 
tures it is difficult to penetrate deeply into this collisional 
hydrodynamic regime as three-body molecule-formation 
will give rise to fast decay of the samples [|. There- 
fore, also the transition region between collisionless and 
hydrodynamic conditions is of substantial practical im- 
portance. 

The hydrodynamic flow of classical fluids was described 
as early as 1755 by the equation of motion of Euler 
The opposite limit of collisionless flow is equally well un- 
derstood since the work of Maxwell and Boltzmann and 
the investigation of rarefied gas dynamics around the 
turn of the last century 0- The transition regime be- 
tween collisionless and hydrodynamic conditions deserves 
special attention as the crossover behavior is often non- 
intuitive as was already noted by Knudsen in 1908 £|. 
With the availability of trapped ultracold gases there is 
a renewed interest in the collisional hydrodynamics. For 
non-degenerate quantum gases in harmonic traps the ab- 
sence of the familiar wall-boundary condition of zero hy- 
drodynamic flow at the sample edges gives rise to a very 
close phenomenological similarity with the superfluid hy- 
drodynamics of Bose-Einstein condensates 0, IE IE 0- 
Collisional hydrodynamics also has to be considered in 
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two-component Fermi gases near inter-component Fes- 
hbach resonances, where the intercomponcnt scattering 
length is tuned to large values in order to optimize ther- 
malization [E 13, E3, Ell 

The onset of collisional hydrodynamics was first ob- 
served at MIT in measurements of the damping and fre- 
quency shifts of the low-lying M — quadrupole shape 
oscillation of cigar-shaped samples of the 23 Na quantum 
gaSjJust above the Bose-Einstein transition temperature 
T c [12]. Similar results were obtained at the ENS-Paris 
with clouds of metastable triplet helium (He*) A 
demonstration of the collisional crossover was given at 
JILA by measuring, for varying density, the damping 
of the center of mass oscillations of two distinguishable 
clouds of 40 K, passing in anti-phase HI EE]. At AMOLF 
we showed how hydrodynamic conditions affect the BEC- 
formation process in elongated samples and can give rise 
to substantial shape oscillations of the condensates being 
formed [1(1 Il7j . Further, hydrodynamic conditions were 
shown to give rise to anisotropic expansion of thermal 
Bose gases after switching off the confining field, which 
has important consequences for time-of-ffight thermom- 
etry [la . Il9l | . Hydrodynamic effects were observed more 
pronouncedly in the expansions of two-component Fermi 
ases tuned near an inter-component Feshbach resonance 
IE El ■ Also the investigation of the macroscopic dy- 
namics of two-component Fermi gases in the BCS-BEC 
transition region requires detailed understanding of the 
hydrodynamics [U El El E| . 

In this paper we study the crossover from collision- 
less to hydrodynamic conditions in non-degenerate clouds 
of 87 Rb by measuring both the frequency shift and the 
damping of the low-lying M = quadrupole shape oscil- 
lation as a function of density. In accordance with theory 
|E IE, , the frequency shifts down from 2u> z in the colli- 
sionless regime to 1.55w z for collisionally hydrodynamic 
clouds, with uj z the axial frequency of our trap. Most of 
the shift occurs over a narrow range of densities around 
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the crossover density, where the mean-free-path becomes 
comparable to the axial size of the sample. At this den- 
sity also the strongest damping is observed. All our data 
were taken for temperatures T > 2T C to avoid precur- 
sor phenomena close to the BEC transition j24|. Hence, 
although the collisions are quantum (i.e. s-wave), the 
gas is statistically classical. As we observed a tempera- 
ture dependence of the oscillation frequencies, we derive 
theoretical expressions to include the influence of trap 
anharmonicities, which cause this effect. These expres- 
sions allow numerical evaluation for regular potentials. 
Further, they allow us to derive convenient analytic ap- 
proximations that apply to any elongated Ioffe-Pritchard 
trap. 



II. BACKGROUND 

For quantum gases well-above the degeneracy temper- 
ature, all oscillatory modes are solutions to the classical 
Boltzmann equation [Hl25| 

^+v-V r / + --V v / = J coll [/], (1) 
at m 

where / = f(t,r,~v) is the phase-space distribution- 
function with r —{rj} — (x,y,z) and v —{vj} the po- 
sition and momentum vectors, m is the atomic mass and 
F(r) = — V r [/(r) the force of the trapping potential U(r); 
Icoi\ [/] is the classical collisional integral. For the case of 
s-wave collisions with energy-independent cross section 
a it takes the form 



/coll[/] 
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d^dSl' |vi - v| [f[ f - h f] (2) 



and describes, for a given position and time, the effect of 
elastic collisions between a pair of atoms with initial ve- 
locities v and vi and final velocities v' and . The solid 
angle fi' gives the direction of the final relative velocity 
with respect to the initial one. 

For isotropic harmonic traps the normal modes are 
multipoles of order (L;M) p. Oscillations in the dipole 
mode (L = 1) are commonly used for measuring trap fre- 
quencies by observing the motion of the center of mass 
of trapped clouds as a function of time. In harmonic 
traps these oscillations are undamped, since for any pair 
of atoms also their center of mass oscillates at the trap 
frequency uj. As noted in ref. Boltzmann obtained 

in f 897 the surprising result that for isotropic harmonic 
traps also the monopole (or 'breathing') mode [L = 0) 
is undamped, oscillating at frequency 2lu, independent of 
the density. The next normal mode solutions are shape 
oscillations (L > 2). In the hydrodynamic limit they are 
(like the dipole mode) both irrotational and divergence- 
free 0, , and therefore also undamped. They oscillate 
at frequency \J~Llo 26] . In the collisionless regime they 
are again undamped but oscillate at frequency L uj. This 
difference in f requ ency results in damping in the tran- 
sition regime |27| due to collisional relaxation towards 
equilibrium. 



For elongated harmonic traps, with axial direction z 
and radial coordinate p — (x 2 + y 2 ) 1 ^ 2 , we distinguish 
three (L = 1;M = 0, ±f) dipole modes, oscillating un- 
coupled and undamped at frequencies uj z and uj p . In the 
hydrodynamic limit, the monopole mode is coupled to 
the (L = 2; M = 0) quadrupole mode. Decoupling in 
terms of irrotational solutions yields 0, IE 



± 



Auji 



25^ + 16c4 - 32uj 2 p uj 2 } 



(3) 



In the experiment described in this paper, we study the 
low-lying M = coupled monopole-quadrupole mode, 
corresponding to the minus sign in Eq. Q. In this mode 
the radial size oscillates in anti-phase with the axial size. 
For shortness we will refer to it as the 'quadrupole' mode 
with frequency loq in all regimes, although in the colli- 
sionless limit the axial and radial motion decouple and 
the overall behavior is to be considered as a superpo- 
sition of 'ID-breathing modes', showing dephasing be- 
havior. This dephasing can be avoided by exciting a 
pure axial oscillation. As follows directly from Eq. J2J, 
in the limit of very elongated clouds {uj p 3> uj z ) we have 

LUQ = y/l2/5iJ z 1=3 1.55 UJ Z . 

The transition regime is less obvious. Describing the 
oscillation phenomenologically by e~ %U3t the crossover 
takes the form [27j 



''hd 



1 — i u>f 



(4) 



where i 



is the complex quadrupole frequency 



for a given thermal relaxation time f; Whd and ll> c \ are 
the (real) frequencies of this mode in the hydrodynamic 
(uj't <C 1) and collisionless (uj't ^S> 1) limits, respectively. 

To have an intuitive picture of the solutions of Eq. Q 
one can separate the real and imaginary parts of uj and 
make the identification (uiq = ui' , T = —u>"). For 
T/luq <C f we can approximate the imaginary part of 
the solution by the convenient form 
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(5) 



which underestimates the maximum damping by 23%. 
The solution for the real part of Eq. @ can be heuristi- 
cally, but fairly accurately (±0.3%) described by 

- ^hd + (w c i - ^hd) (2/7r) arctan(w^f 2 ) . (6) 

We define the 'cross-over point' as the point where max- 
imum damping occurs and the frequency is at the inter- 
mediate value ujq — 1/2 (lo c i + cuhd)- From Eqs. JSJ and 
© this is seen to occur at the relaxation time f = tq, 
where 2lu z tq = 1. 

Eq. Q can be obtained from the Boltzmann equation 
in the relaxation time approximation pel |2jj , where the 
collisional integral is replaced by 



/coll[/] 



(/-/ie)/r. 



(7) 
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Here r is the relaxation time and f\ e = f\ e (t, x, v) the lo- 
cal thermal distribution, which has an isotropic momen- 
tum distribution [2i| • For harmonic traps one has [2^] 



V2 



where r„ 



V2 



n v th a 



(8) 



is the elastic collision rate at the trap center [25|. with 
vth — (SksT/inn) 1 / 2 the mean thermal velocity at tem- 
perature T, no the central density and a the elastic- 
scattering cross section. 

To arrive at Eq. Q the relaxation time has to be renor- 
malized by a factor that depends both on the cloud shape 
and the mode considered. For the M = quadrupole 
mode in very elongated harmonic traps one finds f = 

6/5r mim. 



III. EXPERIMENT 

In our experiment we typically load 10 10 atoms from 
the 87 Rb source described in ref. j^H into a magneto- 
optical trap. After an optical molasses stage we opti- 
cally pump the atoms into the fully stretched (5<Si/2, F — 
2,m,F = 2) hyperfme state and transfer the cloud into 
a Ioffe-Pritchard trap with frequencies u) z /2ir — 7 Hz 
and cjp J 2tt = 8 Hz and central held Bo = 37 Gauss. 
Any remaining population in the mp = 1 magnetic sub- 
level is removed by gravitational sag. Then, we radi- 
ally compress the cloud, changing the trap parameters 
to us p /2ii — 19 Hz at Bq = 8 G. After a thcrmaliza- 
tion time of 100 ms we add, in a linear ramp over 0.5 
ms, a magnetic held B m = 487 mG, rotating at a fre- 
quency of v m = 7 kHz orthogonally to the trap axis, us- 
ing the approach described in ref. [12 ■ This gives rise 
to a Time- Averaged- Potential (TAP) field with offset 

2 , D 2\l/2 



0,m 



[Hi 



BqY'" and frequencies 



(1 + 0.5 6 2 ) 1 / 2 
9 (I + 6 2 ) 3 / 4 ' 
1 

(l + b 2 ) 1 ^' 



(9) 
(10) 



where b = B m /Bo is the modulation depth |32| . 

We continue the compression to Lu z ^ m /2ir — 16.8 Hz, 
u! p>m /2n = 474 Hz and B , m = 634 mG (B = 406 
mG, b = 1.2). Then, we cool the sample by forced rf- 
evaporation to the final temperature of a few microkelvin. 
After reducing the density to the desired level by laser de- 
pletion |33| , the sample is thermalized during plain evap- 
oration periods of up to 2.5 s , which is sufficiently long 
even for our lowest densities. We then raise the rf-shield 
energy by a factor of 7 to avoid evaporation losses during 
the measurements. 

It is important to note here that the harmonic range 
Po,m where Eqs. JHJ and IjlOjl hold is proportional to the 
amplitude of the rotating field: po, m = B m /a, where a 
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FIG. 1: Typical quadrupole oscillation traces for fitting fre- 
quency and damping; raw data acquired at a temperature of 
2 fiK: a) Low density data (no = 1.7x 10 11 cm -3 ) with fit func- 
tion (grey line) resulting in loq/luz = 2.0 and T/ui z = 0.01. 
Each point represents the average of 3 absorption images, 
b) High density data (no = 1.1 x 10 14 cm -3 ) with fit function 
(grey line) resulting in uiq/lu z = 1.6 and T/ui z — 0.22. Each 
point represents the average of 15 phase contrast images. The 
error bars represent the error in the mean. 



is the radial gradient of the Ioffe quadrupole field. For 
regions outside the harmonic range the frequencies revert 
to the unmodulated ones |34| . This implies a minimum 
required value for B m to assure that the harmonic radius 
of the TAP held exceeds the thermal size of the sample. 



A. Excitation of the quadrupole mode 

To excite the quadrupole mode we remove the mod- 
ulation field B m and observe the oscillation in a static 
potential. The advantage of the TAP approach is the 
rapid switching between trap frequencies, which is pos- 
sible because both B m and Bq are generated by trim 
coils. The main currents of the Ioffe-Pritchard trap re- 
main untouched. Further, this approach offers definite 
knowledge of phase and amplitude. After transfer into 
the static potential the cloud starts to oscillate inwards 
as a cosine function with zero phase offset as can be seen 
in Fig.HJ 

As we remove the modulation, we simultaneously in- 
crease the central field to Bq — 900 mG in order to keep 
lu p constant. The procedure is done with a linear ramp 
of duration t sw = 230 us. This is slow enough to avoid 
switch-off depolarization and still much faster than the 
axial oscillation time, uj z <§; r^j <§; ^>Larmor- Thus, 
after switching u) z , the gas finds itself diabatically in an 
axially tighter potential. The axial trap frequency has 
increased to lu z /2tt = 21.1 Hz, which changes the aspect 
ratio to uj p /uj z « 23 and puts us well into the elongated 
trap limit of Eq. © . 

In this way we excite a pure axial oscillation, at least 
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in the collisionless limit. In the hydrodynamic limit, in 
principle both the low-lying and the high-lying monopolc- 
quadrupole modes could be excited. However, since even 
at our highest densities, radially we remain in the col- 
lisionless regime, the high-lying mode cannot be excited 
due to the lack of coupling. 

If the extent of the cloud prepared in the TAP- 
modulated magnetic field reaches significantly beyond 
po,m, its density profile deviates from a Gaussian. Trans- 
ferring that distribution in the described way into the 
static magnetic field leads to excitation of higher modes. 
If additionally the thermal size of the cloud also exceeds 
the harmonic range of the static potential (which is not 
related to po. m ), both anharmonicities will add to the ex- 
citation of higher modes. However, these modes oscillate 
at much higher frequencies and damp accordingly faster 
than the quadrupole mode under investigation. 

For our highest-density samples, together with the con- 
dition T > 2T C , we have po, m ^ 10 A*m, which implies a 
required TAP amplitude B rn > 350 mG. The value of 
B m = 487 mG, used in the experiment, represents our 
technical limit, and corresponds to the harmonic 1/e-size 
of a thermal cloud at a temperature of 9 pK. To assure 
that fitted values for frequency and damping are unaf- 
fected by higher modes, we neglect the first two cycles of 
oscillation traces acquired at temperatures above 7/iK, 
and the first cycle for traces acquired above ApK. For 
lower temperatures also the first cycle is analyzed. Note 
that the precise reproducibility of the starting phase of 
our oscillations allows this procedure without degrading 
the quality of the fits. 



B. Description of the trapping field 

During the observation of the quadrupole oscilla- 
tion the cloud resides in a potential given by U (r) = 
Pb[B(y) — Bq] + mgy for the chosen Zeeman-level in 
this experiment. Here ps is the Bohr magneton and g 
the gravity acceleration along the vertical direction (y- 
direction). For elongated Ioffe-Pritchard traps the mod- 
ulus of the trapping field B(r) is accurately described by 

mm 

B(x,y,z) = \J (B + f3z 2 ) 2 + a 2 (x 2 + y 2 ) + 4a/3xyz, 

(11) 

where Bo = 0.9 G and a = 353 G/cm are defined above 
and 2/3 = 274 G/cm 2 is the axial curvature. To our 
knowledge a and f3 were constant throughout the mea- 
surements to within 0.1%; Bq was monitored to be con- 
stant to within 1%. Expanding Eq. around the trap 
center and keeping the leading non-linearities [3t3 |. the 
potential can be written as 

U(p 1 z) = \m[u 2 z z 2 (\-^p 2 /pl)+ 

+ c> 2 (l-ip 2 /Po)] + -- - , (12) 



where muo 2 — 2u B /3, muj 2 = and 2 / B and po = B /a — 
25 pm the harmonic radius |37j ■ 



C. Detection procedure 

Two imaging methods are used to observe the oscil- 
lations. For our highest-density samples we use phase- 
contrast imaging with red-detuned light. For densities 
tiq > 5 x 10 13 cm~ 3 a proper contrast is obtained at a 
detuning of —3 GHz, where the detection is essentially 
non-destructive |2^. This allows us to register the os- 
cillations in a sequence of 31 images at 5 ms intervals, 
taking adva ntag e of the fast 'kinetics' imaging mode of 
our camera |39j. For lower densities the phase contrast 
method cannot be used because, at the (smaller) detun- 
ings required to maintain adequate phase contrast, pho- 
toassociation losses disturb the measurements 4Q] . 

For densities no < 5 x 10 13 cm~ 3 , we used repetitive 
absorption imaging on the (5Si/ 2 , F = 2) <-> (5P 3 / 2 , F — 
3) transition (Z^-hne) (4^, varying the holding time of 
the cloud after excitation of the oscillation. The images 
were taken in situ, just before releasing the cloud from 
the trap 01 • We apply the usual method of background 
subtraction and level- normalization to process the images 
m>mi- To retrieve the column density profile n 2 (y,z) 
and the axial and radial Gaussian 1/e-sizes L e and R e , 
we fit a 2-dimcnsional Gaussian expression to the optical 
thickness distribution of our images. The central density 
follows with no = n 2 (0, 0)/ \JitR 2 and, with Eq. JHJ, the 
relaxation time can be expressed as 

" 2T= 2^^0y- (13) 

Note that this expression does not depend explicitly on 
the gas temperature. The collision cross section is a = 
87ra 2 in the zero temperature limit and is calculated with 
the value a = 98.98(4)ao for the s-wave scattering length 

m. 

To acquire sufficient statistics, at least 30 images are 
taken to retrieve one oscillation trace for a given density 
and each trace is acquired at least 3 times. Because the 
crossover happens over a narrow range of densities, great 
care was taken to reproduce the initial conditions from 
shot to shot. This is done by adjusting the density using 
laser depletion in a feedback loop with the experimental 
result of the previous shot m| . Although this procedure 
increases the shot to shot fluctuations, long-term drift is 
virtually eliminated. With this procedure the atom num- 
ber could be long-term stabilized within a standard devi- 
ation of better than 1%. By fitting the expression for an 
exponentially damped cosine function to the trace (see 
Fig-Di we retrieve the experimental values for the fre- 
quency u>q and damping rate T of the quadrupole mode 
for the selected density. 
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D. Accuracy of density and temperature 
determination 



The absolute accuracy of rt2(0,0) is estimated to be 
~ 30% The phase contrast images are calibrated 

against absorption images of expanded clouds taken 15 
ms after release from the trap at zero detuning. This pro- 
cedure presumes the conservation of atom number during 
the expansion. 

In our analysis we account to leading order for the cor- 
rections associated with trap anharmonicities. For tem- 
peratures much lower than the harmonic temperature 
To = Ms^o/fes = 60/xK, Eq. (JT2J becomes sufficiently 
accurate to describe the cloud shape. In this limit the 
column density on the trap axis (to leading order in the 
x-integration) can be expressed for z 2 -C 2k B T/moj 2 as 



n 2 (0, z) ~ n 2 {0, 0) exp 



2k B T 



.(1--T/T„) , (14) 



where T/T = (x 2 ) /pi = k B T/fi B B with (a; 2 ) 
k B T/mui 2 the variance of the thermal distribution of the 
cloud along the x-axis in the harmonic limit. 

From Ea. (|14|l we estimate the 1/e-axial-size L e that 
will be obtained by fitting a Gaussian to the axial column 
density profile of the cloud, L e = L/(l — jT/Tq) with 
L defined by L 2 — 2k B T /muj 2 49]. The temperature 
follows with the expression 



k B T 



±mtu 2 z L 2 e (l-±T/T ). 



(15) 



Hence, for a temperature of 6 /iK the harmonic approx- 
imation overestimates the temperature by ~ 5%. The 
correction in the central column density is smaller. Nu- 
merically we established that the fit of a 2D Gaussian 
underestimates the central column density by ~ 1.4% at 
T/Tq = 0.1. As these corrections are small, there is no 
need to go beyond the leading order of anharmonic cor- 
rection to retrieve these quantities. For measuring oscil- 
lation frequencies the situation is different because these 
can be measured to high precision. 

Mean field broadening of the distribution is small [5(| • 
Calculating the variance (z 2 ) = ^L 2 using the recur- 
sive expression for the density to first order in mean field 
U mi (r) = 2u n(r), leads for T < T to 



1 



2 T i 



—muj^L 
2 



k B T + E, 



mi 1 



(16) 



where E m{ = Vq J n 2 (r)dr/ J n(r)dr = v n /\/8 is the 
trap averaged interaction energy with vq = (AttH 2 /m)a 
the interaction coupling constant 1] . Equivalently, treat- 
ing the mean field as an effective potential we may write 



1 



,2r2 



(17) 



where £ = E m {/ (k B T + E m {) is the mean field correction 
constant. For the data point with the highest mean field 
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FIG. 2: Scaled frequencies Cjo/^Jd for the axial dipole mode 
versus temperature. The grey line corresponds to the evalu- 
ation of Eq. (1241 . The leading slope (dashed line) is given 
by Eq. 12511 . Open circles: phase contrast measurements; 
closed circles: absorption imaging measurements, both ac- 
quired with thermal clouds. Diamonds: measurements with 
Bose-Einstein condensates, using absorption imaging. Open 
squares: results acquired at ENS-Paris |5lj |. 



= (n = 1.1 x 10 14 cm~ 3 , T = 2 uK) we calculate £ = 0.007. 



Therefore, mean field corrections are at least one order of 
magnitude smaller that the anharmonic corrections and 
are discarded in this paper. 



IV. ANHARMONIC FREQUENCY SHIFTS 

As we operate at temperatures well above T c , we pay 
special attention to the issue of trap anharmonicities. We 
follow the path of argumentation as presented in [2^, 
to derive expressions for the anharmonic shifts. These 
are both temperature and mode dependent and can also 
depend on the density. The expressions are suitable for 
numerical evaluation provided the first and second spatial 
derivatives of the trapping potential are known. 

To describe the dynamical evolution of an observable 
X = x( r i v ) if is multiplied by Eq. (JIJ and averaged over 
the phase space 



dt 
where 



VrX) 



(F • V vX ) 



(x) - (x)ic 



(x)(t) 



1 

N 



x(r, v)/(t, r, v)d 3 rd 3 v, 



(18) 



(19) 



with N the number of atoms. By choosing the correct 
set of observables, it is possible to obtain a closed set of 
equations that describes the dynamics of these observ- 
ables. 



A. Dipole mode (L = 1) 

To investigate the effect of trap anharmonicities on the 
dipole mode oscillation, we make the Ansatz 



/(*) r , v ) = fo(n - a i} m - hi), 



(20) 
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where / (r, v) = C exp(- (mv 2 /2 + U(rj) fk B T) is the 
equilibrium distribution function with C the normaliza- 
tion factor and = fli(i). We choose x — v i an d obtain 
the following set of equations 



Tit 



1 



(v^ (fi(r))=0, where 



and (F(r)) = (F(r + a)) c 



(21) 



(22) 



Analogously to Eq. (|19|l we denote with (x)o the aver- 
age on the phase space using the equilibrium distribution 
/o(r, v). Expanding up to first order around the equilib- 
rium position a.; = 0, we obtain 



1 

m 



0, 



(23) 
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FIG. 3: Scaled frequencies uoq/ujq for the quadrupole mode 
versus temperature. The grey sector corresponds to Eq. I|M7[I 
evaluated for the cross-over regime (see Sec llVB"3t with the 
collisionless and hydrodynamic limits indicated. The dashed 
lines show the leading slopes given by Eqs. (gOJ and (HJ. All 
data points correspond to collisionless conditions. The dotted 
line is used to scale all quadrupole frequency data to zero 
temperature. 



where U'A = 



Restricting ourselves to potentials 

with (U"j)o = for i ^ j, we obtain for the effective 
frequencies of the dipole modes 



Co 2 



L< - —(Uii)0- 

m 



(24) 



Substituting Eq. ljT2l for the potential into Eq. (|2^|) we 

obtain for the leading anharmonic shift in the z-direction 



&zD ^ - -T/Tq). 



(25) 



This expression is shown as the dashed line in Fig. [2 
The integral in Eq. (|24(l is readily evaluated numerically 
using Eq. 1)11(1 and requires as input parameters only the 
values for a, (3 and B . The resulting curve is shown as 
the solid line in Fig.|2 The curve follows the trend of 
our measurements of center-of-mass oscillations as well 
as data obtained in Paris [5l| . 

The zero temperature limit of uj z d is largely fixed by 
measurements with Bose-Einstein condensates, which re- 
produced within 1% over a period of one year. Its value 
is used to calibrate oj z and the related /^-coefficient. We 
have no explanation for the remaining deviations for 
the points taken with thermal samples at higher tem- 
peratures [52^|. We cannot trace them back to insuffi- 
cient mechanical or electronic stability of our trap. Non- 
exponential contributions to the damping may account 
for a systematic error in the frequency, but should be less 
than 1%. We speculate that possibly the temperature de- 
terminations of the phase contrast measurements could 
be affected by a molecular contribution to the phase con- 
trast, which tends to narrow-down the distribution and 
results in an underestimated value for the temperature. 
This results from the distribution of pairs, that can pho- 
toassociate, which is proportional to the square of the 
atomic density. 



B. Surface modes (L — 0, L = 2) 

In order to calculate the anharmonic shifts for the 
breathing and the two quadrupole modes we make the 
Ansatz 



/(*,r,v) 



1 



/o U? ' (26) 



where hi and 9i are time-dependent variables. The 
parameters bi take into account shape deformation of 
the density cloud whereas the parameters 9i allow an 
anisotropic momentum distribution which is crucial to 
calculate the correct frequencies. We choose x = v i r i 
and obtain the following set of equations 

k{r*) - |(« 4 2 )o - ^(FifarAnh = 0. (27) 

We impose the stationary solution and find the relation 

1 



K)o = (Fir^o. 

m 

Then choosing \ = ( v i — (bi/bi)ri) 2 yields 
x Bi-B 



(28) 



(29) 



where 0= (E* W- 

Let us now focus our attention on two extreme regimes. 



1. Collisionless limit 

In the collisionless limit (r — > oo) we obtain the rela- 
tion 



1 

w 



(30) 
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and finally, 

b,- 



1 (i ? i(r 7 -)ri)o {Fi{b 3 r 3 ) ri ) _ Q 



6f m(r?) 



m(rf) 



Linearizing these equations and looking for solution of 
the form e~ lult we obtain the three frequencies. In order 
to do this explicitly we define the quantities 



2. Hydrodynamic limit 

In the hydrodynamic regime (r — > 0) the local equilib- 
rium is always maintained, which implies that 9i = 9. In 
this case we obtain the relation 



1 



At 



3 (nUDo 

m (r 2 ) 



2/3 ' 



where U[ = Aij — for i ^ j and 



Bi 



1 (nrjUtyo 
m (r?) 



(32) 



(33) 



and therefore 



(Fi(r-,-)r z ) (Fjjbjrjy^o 



MILA) 2/3 m(r?) 



m(r?) 



(41) 



= 0. (42) 



where I7| • - r= 

have to solve 



; note that, in general, B^ £ B 3i . We fine the A ma trix as 



By linearizing around the equilibrium we find the fre- 
quencies for the M — ±2 modes and the two M = 
monopole-quadrupole modes. In this case we have to de- 



\A + B-Cu 2 I\ = 0, 



(34) 



where I is the identity-matrix, in order to obtain the 
frequencies. For the quadrupole modes with M = ±2 we 
find 



J Q2 



AxX A X y -\- B XX B X y 



whereas for the modes with M = we have 



<D B 2 = (i + \J [i 2 - A 

- 2 



= /x - v 7 ^ 2 - A, 



(35) 

(36) 
(37) 



where 



M — (A X x ~\~ A X y + v4 zz + + -B^y + £U/2 (38) 

A = (Ajjj; + A X y + -B^j; + B X y)(A ZZ + ) 



(39) 



Here we used (JJ"j)o = for i ^ j. An analytic approx- 
imation for the leading anharmonic shift of the M = 
quadrupole mode is obtained by substituting Eq. (|12|) for 
the trap potential, 



-T/T ), 



(40) 



where wq is the frequency in the harmonic limit. Note 
that at this level of approximation the relative shift coin- 
cides with that of the dipole mode. The result of the nu- 
merical averages based on Eq. (|ll|l is shown as the lower 
solid line in Fig-EI The lower dashed line corresponds to 
the leading shift given by Eq. (|40p . 

Comparison with the experimental points in Fig-El 
shows agreement as far as the trend of the shift is con- 
cerned but a systematic deviation for the slope. This dis- 
crepancy can be eliminated by presuming that our axial 
trap frequency lo z is underestimated by 1.5%. However, 
such a correction cannot be justified on the basis of the 
limited set of data for the dipole mode |5^ . 



An — 



5 (nUDo 

3m (r 2 )n 



(43) 



Note that A^ does not depend on j. The matrix Bij is 
the same as in the collisionless case. Solving the deter- 
minant Eq. (JHU leads again to Eqs. J3U), (j2EJ) and IpTfjl 
for the frequencies and Eqs. (|38(l and (|39|l for fx and A. 
Only the expressions for the matrix elements A^ have 
changed. Substituting Eq. (|12(l for the trap potential we 
find for the leading anharmonic shift of the hydrodynamic 
M = quadrupole mode 



u Q ~ uj q (\ - -T/T ), 



(44) 



which has a slightly weaker slope than in the collision- 
less case. The result of the numerical averages based on 
Eq. (jl 1|> are shown as the upper solid line in Fig-EJ The 
upper dashed line corresponds to the leading shift given 
by Eq. gU). 

A comparison with experiment requires densities hq > 
4x 10 14 cm -3 at a temperature T = 4 /iK, to have 2ui z f < 
0.1, which is about three times our maximum density. At 
our highest density of no = 1.3 x 10 14 cm~ 3 , we calculate 
a 3-body decay rate of N/N = 2\/3 Ln 2 , « 1 s _1 , with 
L = 1.8(5) x 10 _29 cm 6 s _1 the three-body rate constant 
in the Bose-condensed state . At a three times higher 
density, the decay rate renders the acquisition of data at 
approximately constant density impossible for 87 Rb. 



3. Crossover regime 

In the cross-over region, the same approach can be 
used, but after linearizing, one has to look for solutions 
of the form e~ luJt with a complex u>. For the M = 
modes this leads to the equation 



C[oj] - ~D[w] j ( E[u] - -F[lo] ) = 0, ( 4r, > 
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where C[lu] — lo(lo 2 — lo 2 ib )(lo 2 — w^.q), D[lo] — (lo 2 — 
^L-.b)(^ 2 ~ w hd: Q )> E[u\ = uj(uj 2 - uj 2 1:Q2 ) and F[u] = 
(lo 2 — <^hd-Q2)- Each term represents two equations since 
they contain real and imaginary parts. For an elongated 
cigar-shape trap it is possible to write the frequencies 
in the form of Eq. J2J with rescaled relaxation time f = 
( <jJ ci-B/ a; hd B) r ' wn ich reaches the value f = 6/5 r in the 
harmonic limit. The numerically calculated results of 
temperature induced shift, based on Eq. (fTTfl in the cross- 
over regime is represented by the grey sector in Fig.|3| 

A comparison with experiment is beyond the scope of 
this paper because, after scaling to loq, the two limiting 
cases are spaced by only ~ 1%. Therefore, not only loq 
has to be determined to an accuracy much better than 
1%, but also the scaling parameter loq. In the limiting 
cases the latter is fully determined by the trap frequency 
(i.e. uj c i = 2lo z and lo^ w 1.55w z ). However, in the 
crossover region knowledge of f to much better than 1% 
is required to calculate loq from Eq. © to adequate pre- 
cision. 



V. RESULTS AND DISCUSSION 

We took all our data with the same trap parameters 
and the same excitation procedure, but at various tem- 
peratures. Starting the evaporation with a large atom 
number and using 'tight' trapping parameters we could 
reach high densities and thus study the full crossover. 
However, this choice for a 'tight' trap made us sensi- 
tive for anharmonic shifts as discussed in section llVl 54] . 
Therefore, we extrapolate all frequency data to the zero- 
temperature limit (Cjq — ► loq) using the dotted curve in 
Fig-El This yields the best estimate for the value in the 
harmonic limit of our potential. The correction curve is 
based on the temperature dependence observed for our 
data in the collisionless regime (2lo z t > 10, see Fig.01, 
where we may presume loq — 2lo z . In this way system- 
atic deviations of our results from the curves in Fig. 01 
and Fig.EJ} were substantially reduced. 

In Fig.0] we plot the observed, scaled damping rates 
T/lo z versus the extrapolated quadrupole frequencies nor- 
malized to the axial trap frequency, loq/lo z . The drawn 
curve corresponds to the crossover expression, Eq. flljl 
with cj c i = 2lo z and Whd = 1.55w z . Plots of the same 
experimental data and the exact solutions of Eq. Q sep- 
arately against 2lo z t are given in Fig.[SJi for T/lo z and 
Fig.|St> for uq/lo z . 

From the damping results (Fig.[S^,) we obtain 2lo z tq = 
1.0(1) for the experimental value of the crossover point. 
Given the 30% absolute accuracy of our density determi- 
nation this agreement is fortuitously good (see Sec lTT|) . 

The determination of the crossover point from the fre- 
quency crossover behavior is less straightforward, be- 
cause errors in the temperature determination add to 
the error in the extrapolated frequency loq. Further, 
as the frequency corrections are all positive, they af- 
fect the determination of 2lo z ?q from Fig.JSJx For the 



crossover region (0.1 < 2lo z t < 10) the average applied 
frequency correction was Sloq/lo z = 2.6%, which changes 
the experimental value for 2lo z tq by a factor 0.85 to yield 
2lo z tq — 1.0(1) j5f|. This value coincides with the re- 
sult obtained from the damping data and shows that our 
results are self-consistent. 

In the collisionless regime anharmonicitics can give rise 
to dephasing induced damping. These effects were not 
corrected for as they do not affect to leading order the 
determination of the crossover point. Here we briefly 
comment on these effects. Roughly, one may argue that 
for a given anharmonic spread Sloq in frequencies the 
dephasing time St will be given by 5loq51 w 2ir. Hence, 
the dephasing related damping rate is ST = 2n/dt w 
Sloq. The cluster of data points at 2lo z t 30 in Fig.|5^ 
best illustrates the significance of the correction as they 
were taken at the highest temperature (9 /xK). For these 
points the anharmonic frequency shift Sloq is ~ 3.5% (see 
Fig.0. With 8T/lo z ~ 5lo q /lo z = 0.035 this suggests 
that the anomalously high damping rates observed for 
these data points (Fig.lSJi) may be entirely attributed to 
dephasing effects. Near the crossover point the collisional 
damping is much faster and dephasing corrections may 
be neglected ST/T ~ (1/2) (Sloq/T) 2 « 10~ 3 . 

We also verified that our shot-to-shot variations in 
the density have a negligible effect on the measured 
damping rate. The frequency shifts fastest at the den- 
sity of the crossover point, where S(loq/lo z ) /5(lo c \t) = 
(2/ir) (u> c \ — LOhd) /w z ~ 0.3 as follows directly by tak- 
ing the first derivative of Eq. JSJ) with respect to f . As 
f scales inversely proportional to the central density, a 
1% variation in atom number results (at constant tem- 
perature) in a ~ 0.3% variation of the frequency, which 
is much smaller than to the one considered above and 
therefore also negligible. 




FIG. 4: Damping (raw data) versus frequency (extrapolated) 
for the quadrupole mode. Open circles: data acquired at 
MIT 0, ENS [3 and AMOLF 0. The solid line is cor- 
responds to the crossover expression Eq. @. The data left 
of the straight line are obtained with Phase Contrast (Ph.C.) 
imaging, those on the right with Absorption (Abs.) imaging. 
The error bars represent a 95% confidence interval of the fit. 
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FIG. 5: Damping (raw data) (a) and frequency (extrapolated 
to harmonic limit) (b) of the quadrupole mode versus scaled 
relaxation time. The solid line in figure a) corresponds to 
the imaginary part (F = — uj") of the solution of Eq. @; in 
figure b) to the real part (ujq = uj'). The crosshair marks the 
location of the crossover point. The vertical error bars are 
identical to Fig.|lJ the horizontal ones represent the standard 
deviations. 



VI. SUMMARY AND CONCLUSIONS 

With Fig.0] and Fig.[5] we obtain good agreement be- 
tween experiment and the crossover theory. The fre- 
quency shifts down from 2uj z in the collisionless regime to 
1.55w z for collisionally hydrodynamic clouds, with ui z the 
axial frequency of our trap. Most of the shift occurs over 
a narrow range of densities around the crossover point. 
The damping rate peaks over the same range of densi- 
ties. The determinations of the crossover point from the 
frequency and the damping behavior agree within 10%, 
2w z fo = 1.0(1). The agreement with theory is limited 
by a 30% absolute uncertainty in density. Further, we 
present a theory and experimental evidence for anhar- 
monic frequency shifts. The theory allows numerical eval- 
uation for potentials with known first and second spatial 
derivatives. We show that for elongated Ioffe-Pritchard 
traps knowledge of the central field Bo suffices to calcu- 
late the leading anharmonic shifts with simple analytic 
expressions. 
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